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Abstract 

A hyperbolic BC{n) Sutherland model involving three independent coupling constants 
that characterize the interactions of two types of particles moving on the half-line is 
derived by Hamiltonian reduction of the free geodesic motion on the group SU{n, n). The 
symmetry group underlying the reduction is provided by the direct product of the fixed 
point subgroups of two commuting involutions of SU{n,n). The derivation implies the 
integrability of the model and yields a simple algorithm for constructing its solutions. 
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1 Introduction 



The integrable many-body systems discovered by Calogero and Sutherland around 40 years 
ago still enjoy extraordinary popularity due to the wealth of their physical applications and 
connections to mathematics, which are described in the surveys |8l |T2l |15l [17]. In correspon- 
dence to the many variants of these systems (associated with different interaction potentials, 
root systems, relativistic deformations, inclusion of two types of particles, and so on) there ex- 
ist also several approaches to studying them. The systems based on trigonometric/hyperbolic 
interaction potentials are usually called Sutherland type and here we study a particular case of 
such systems in the classical Hamiltonian reduction framework, reviewed for example in |12] . 

A Sutherland type integrable system describing the interaction of m 'positively charged' and 
(n — m) 'negatively charged' particles was first introduced by Calogero by means of shifting 
the positions of m out of the n particles by i^. This trick converts the repulsive interaction 
potential sinh~^(gj — q^) into the attractive potential — cosh~^(gj — q^) between the particles 
of opposite charge (indexed say by 1 < j < m < < n). Then Olshanetsky and Rogov [lOj 
derived the Calogero-Sutherland model by Hamiltonian reduction of free motion on an affine 
symmetric space. The dynamics of the model and its relativistic deformation was analyzed in 
detail by Ruijsenaars [T4], motivated mainly by the relation of this problem to the interaction 
of sine-Gordon solitons and anti-solitons. 

In a little noticed paper Hashizume [7] generalized the Olshanetsky-Rogov derivation and 
thereby proved the integrability of a family of hyperbolic Sutherland models associated to the 
so-called root systems with signature. In the case of the BC{n) root system, his model involves 
two types of particles moving on the half-line with the interaction governed by two independent 
coupling parameters. 

It is well-known that integrable BC (n) Sutherland models involve in general three arbitrary 
couplings corresponding to the three different root lengths. However, it has been explained 
only rather recently, by Pusztai and one of us [S], how the three couplings arise in the setting 
of Hamiltonian reduction. In the present paper, we generalize the result of [S] and derive the 
following BC{n) type Sutherland Hamiltonian, 

1 " 2 2 



2 ^ . ^ cosh [qj — qh) cosh (o,- + qi.) 

J = l l<]<m<k<n 



ij^<ra sinh^(gi - qk) sinh2(g^. + g^) smh^^q^ - qk) smh?{qj + qk) ' 

1 ^ [x — y) 1 ^ xy 1 ^ xy ,^ ^. 

+ 2 ^ sinh2(2g,) + 2 ^ ^^^^) ~ 2 ^^^h^' ^ ' ^ 

where «; > 0, a; and y are real coupling constants. If (x^ — y"^) ^ 0, then energy conservation 
ensures that the corresponding dynamics can be consistently restricted to the domain where 

gi > g2 > ••• > > and qm+i > qm+2 > ■■■ > > 0. (1.2) 

Supposing also that xy > 0, the Hamiltonian (11.11) describes attractive-repulsive interactions 
between m 'positively charged' and (n — m) 'negatively charged' particles influenced also by 
their mirror images and a positive charge fixed at the origin. 
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We shall derive the model (11 .ip by reduction of the free geodesic motion on the group 
G = SU (n, n). Our reduction relies on a symmetry group of the form x G~^, where G+ < G 
is a maximal compact subgroup and G^ < G is the (non-compact) fixed point subgroup of 
a G-involution that commutes with the Cartan involution fixing G+. Hashizume [7j reduced 
the geodesic motion on affine symmetric spaces such as G/G~^, which itself is the reduction of 
the free motion on G associated with the zero value of the moment map of the G~^ symmetry. 
He obtained the model (II. ip with two independent couplings, while we obtain it with three 
arbitrary couplings since we use non-trivial one-point coadjoint orbits of G^ having a free 
parameter (corresponding to y in (II. ip ) to define our reduction. In [5| (see also [H Ej) the 
symmetry group (?+ x G'_|_ was used in an analogous manner to describe the m = case. 

To be more precise regarding the comparison with Ref. [7j, note that the model (II. ip with 
y = was obtained in by using G = U{n, n) and the model with a certain non-linear relation 
between the three couplings was obtained by using G = U{n + l,n). The possible alternative 
reduction treatments of the model (11.11) are briefly discussed also in the concluding section. 

Our derivation implies the Liouville integrability of the model (11.11) in the general case, and 
it also gives rise to a simple linear-algebraic algorithm for constructing the solutions. It could 
be interesting to analyze the dynamics of the model in the future by utilizing this algorithm, 
and to possibly relate it to special solutions in a field theory on the half-line. Further comments 
on open problems are offered at the end of the paper. 



2 Group theoretic preliminaries 

We here fix our notations and recall some group theoretic results that will be needed later. 
To begin, we choose some integers 

l<m<n (2.1) 

and define the matrices 



Qn,r, 



1„ 
1„ 



e gl{2n, C), Im ■■= diag(l^m, -In-m) e gl{n, C), (2.2) 



where 1„ denotes the n x n unit matrix. We also introduce 

Dm := diag(/m, Im) = diag(l„, -ln~m, Im, -In-m) ^ gli2n, C). (2.3) 

We adopt the convention in which the group G := SU{n,n) and its Lie algebra Q := su{n,n) 
are given by 

SUin, n) = {ge SL{2n, C) | g^Q^^ng = Qn,n} (2.4) 

and 

su{n, n) = {V e sl{2n, C) | V^Qn,n + Qn.nV = 0}. (2.5) 
In the obvious n x n block notation the elements V G su{n, n) have the form 



V 



X Y 

Z -xt 



= -Y, = -Z, Q{tT{X)) = 0. (2.6) 
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We consider the commuting involutions of SU (n, n) provided by the Cartan involution 
and the involution F: 



(2.7) 



The fixed point subgroup of © is the maximal compact subgroup G+ < G and the (non- 
compact) fixed point subgroup of F is denoted by G"*". Let 9 and 7 be the corresponding 
involutions of ^ = su{n, n). Using the n x n block notation, the Lie algebra of G+ reads 



I ^ ^ : Xt = -X, Ft = tr(X) = o| , 



(2.8) 



and is isomorphic to s{u{n) © u{n)) according to 



s{u{n) © u{n)) 3 



a 
/3 



^^|^{a,p) := - 



a + /3 a — /3 
a — (3 a + (3 



Correspondingly, the group G+ is isomorphic to S{U{n) x U{n)) via the formula 



S{U{n) X ?7(n)) 9 
which can be written also as 



a 
h 



^ ^(a,6) := - 



a + b a — b 
a — b a + b 



eG, 



g{a, b) = is:(diag(a, b))K-^ with K :-- 



The Lie algebra of G^ is furnished by 



V2 



In In 
In In 



(2.9) 



(2.10) 



(2.11) 



Y 



: Xt = -ImXIm, Ft = tr(X) = 



(2.12) 



and is isomorphic to s{u{m, n — m) ® u{m, n — m)) via the map 



s{u{m, n — m)® u{m, n — m)) 3 



a 

p 



a + /S (a - /S)/™ 

(a — 13) Im{ct + /3)/r7 



(2.13) 

In the above formula u{m, n — m) is realized as the Lie algebra of the nxn matrices satisfying 
the relation 

aU^ + /^q; = (2.14) 

and it holds that 



x(q;,/5) = A'(diag(a,/3))A:-^ with K :-- 



1 

72 



In In 



(2.15) 



Similarly, is isomorphic to S{U (m, n — m) x U (m, n — m)) by means of conjugation by if. 
The eigensubspaces ^_ of ^ and ^~ of 7 having eigenvalue —1 can be displayed as 



X Y 
-Y -X 



: Xt = X, Ft = -f| 



(2.16) 



and 



X 



Y 



X"^ = I XI = -Y 



(2.17) 



Introducing Ql := fl Q'" for any signs s,r E {±}, we can decompose Q as the direct sum of 
disjoint subspaces, 

g = Gz®gi®g^®gi, (2.18) 

which are pairwise perpendicular to each other with respect to the invariant scalar product on 
g defined by 

(2.19) 



{V, w) := -tT{vw), VI/, w eg. 



We have 

gt 



and the appropriate restriction of the map (12.91) gives rise to an isomorphism 

g^ ~ s{{u{m) © u{n — m)) © {u{m) © u{n — m))). 
We shall also use the explicit form of the other subspaces: 

-X^ Y = -I YI = -F^ 

X IfY^X X , Y IfjiY Ifji Y 

■ X = I^XI^ = X\ Y = ImYIm = —Y^ 



(2.21) 



gi 
g- 



X Y 
Y X 

X Y 
-Y -X 

X Y 
-Y -X 



. Y 1 Y T 



(2.22) 



Next, for our purpose we choose a maximal Abelian subspace A of gz, i.e., an Abelian 
subalgebra of g which lies in gz and is not properly contained in any Abelian subalgebra of 
the same kind. It is known [16] that any two choices are equivalent by the conjugation action 
of Gt on gz, and concretely we choose 



A:-- 





' q " 


{" 


-q 



: q = diag(gi, ...,g„), Qk G 



(2.23) 



One can verify that the centralizer of ^ in ^ is given by the direct sum 
C = A®M, M = ld:=i 



d 
d 



: d = diag(c/i, d„), 4 e M, tr(rf) = o| < g+. 



(2.24) 

Denote by A and M the connected subgroups of G corresponding to the Abelian subalgebras 
A and M, respectively. In fact, M is precisely the subgroup of Gt whose elements g satisfy 



gqg 



q, Vg G A. 



(2.25) 
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Furthermore, we call an element q & A regular if the elements g G Gl^_ satisfying the relation 
919"^ = 1 all belong to M. It is not difficult to check that g G ^ is regular in this sense if and 
only if the following conditions hold: 

gj 7^ i = 1, n, (gj — qk){(lj + gfe) 7^ 1 < j < k < m and m < j < k < n. (2.26) 

Choose a connected component Ac of the open set formed by regular elements of A, and 
denote the closure of this 'open Weyl chamber' by Ac- According to general results p^, every 
element g E G can be decomposed in the form 

(7 = (7+eV with qeAc,g+eG+,g+eG+. (2.27) 

The constituent q that enters this decomposition is unique, and if q is regular then the ambiguity 
of the pair {g+,g^) is exhausted by the replacement 

((?+, g+) ^ (g+fi, ii-'g+) G M. (2.28) 

In the the generalized Cartan decomposition (12.271) the open Weyl chamber Ac can be taken 
to consist of the elements q in (12.231) that are subject to the condition 

qi> q2> ••• > gm > and g^+i > qm+2 > ••• > qn > 0. 



Both and possess one-dimensional centres. The centre of is generated by 

Ir 

Ir 



c 



1 







(2.29) 



(2.30) 



and the centre of is spanned by 



In. 

Im. 



(2.31) 
(2.32) 



These elements enjoy the property 

G^eM^f^gX A = /,r. 

The decomposition (I2.27P and the property (12.321) will be important for us in Section 3. 

By means of the invariant scalar product (I2.19p . we can regard Q, and as their 
own dual spaces, respectively. This then also identifies the respective coadjoint actions with 
the adjoint actions. In the next section we shall utilize particular coadjoint orbits of To 
describe them, for any non-zero column vector m G C" define the matrices 



Xiu) 



l[UU 



t 



and ^{u) 



X{u) 
Xiu) 



X{u) 
Xiu) 



(2.33) 



1 

n ' - ■ 2 

Fixing arbitrary real constants n > and x 7^ 0, it is easy to see (cf. (l2.9p -f l^rTT]) ) that the set 

:= {xG^ + ^{u)\ue C", u^u = Iku} (2.34) 

is a coadjoint orbit of G+ of minimal non-zero dimension. The action of gia^ b) G on Ok,x 
takes the form 

g{a, b){xG^ + ^iu))g{a, b)-^ = {xG^ + ^(au)). (2.35) 

Since ^(m) determines u up to an overall U{1) phase, the orbit O^.x can be identified with the 
the complex projective space CP„_i. 

We remark that, for any real constants x and y, xG^ and yG^ represent one-point coadjoint 
orbits of and G+, respectively. 
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3 Hamiltonian reduction 



We shall reduce the free geodesic motion on the group G = SU {n, n) formulated as a Hamil- 
tonian system on the cotangent bundle T*G. We find it convenient to analyze the reduction 
by using the so-called shifting trick of symplectic geometry, which amounts to extending the 
phase space by a coadjoint orbit before reduction [TT]. Specifically, trivializing T*G by right- 
translations and identifying Q* with Q by means of the invariant scalar product, we consider 
the phase space 

P := T*G X a,,. ~ (G X e;) X a,x = {{g,J,0\geG,J eG,Ce a,.}- (3.1) 
The symplectic form on P is given by 

n = nT*G + ^o.,^. (3.2) 

where VLt*g can be written explicitly as 

^lT'G = d{J,dgg-') (3.3) 

while the explicit form of the Kirillov-Kostant-Souriau symplectic form x of the coadjoint 
orbit Ok .a; fl2.34l) will not be needed. The phase space P carries the commuting family of 
Hamiltonians provided by 

Hk{g, J, := ^t^^''), = 1, 2, n, (3.4) 

the first member of which is responsible for the geodesic motion. These Hamiltonians are 
explicitly integrable; the flow of if^ with initial value ((^o, Jqi Co) is readily verified to be 

{g{t), Jit), at)) = (e'^'^g,, Jo, Co) with Vk := J^"'' - ^tr( j2'=-^)l2n. (3.5) 

Note that Hk is real since J^'' satisfies (J^'')"'' = Qn,nJ'^''{Qn,n)~^ , and Vk in (13. 5p belongs to 
Q = su{n, n). 

We introduce an action of the group G+ x G^ on P by sending the pair (r^, h) G (7+ x G^ 
to the symplectomorphism '^ri,h of P operating as follows: 

^,,^((7, J, C) := {'n9h~\ VJV~\ VCV~')- (3.6) 

The Hamiltonians Hk (13. 4p are invariant under this group action, which is generated by the 
equivariant moment map 

$ = ($+,$+) (6;+, 6;+), (3.7) 

$+(^?, J, C) = 7r+( J) + C, <^^{g, J, C) = -n+{g-'jg), (3.8) 

where the projection 7r+ : ^ — t- ^+ is given by means of the decomposition ^ = ^+ © and 
7r+ : ^ ^ ^+ by ^ = g+®g~. 

We are interested in the reduction defined by imposing the moment map constraint 

$ = i/ with u := (O^-yC), (3.9) 
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where ?/ 7^ is a real constant and we refer to (13.71) . The action of the symmetry group x 
preserves the 'constraint surface' 

■= c P. (3.10) 
We require that the constants x and y verify 

ix'-y')^0. (3.11) 

Then the corresponding space of orbits, 

P,ed := Pc/{G+ X G^), (3.12) 

will turn out to be a smooth manifold. According to the general theory [TT], P^d inherits the 
symplectic form fij-ed and the reduced Hamiltonians Hl'^'^ defined by the formulas 

7r*f2i.cd = ^|pc, ^"Hl'"^ = Hk\p^, (3.13) 

where tt : ^ P^ed is the natural projection and ^\p^ is the restriction of Q (13. 2p on P^ C P. 

Remark 3.1. In this technical remark we explain why the space of orbits (I3.12p is a smooth 
manifold. First, we note that the action (13. 6p of x G^ is proper. By definition [11], this 
means that for any sequences {rjn, hn) in G^ x and ((?„, Jn,Cn) in P (with n G N) for which 
{.9ni JnXn) and {q„,h„){.9ni JnXn) ar© both convergent, there exists a convergent subsequence 
of the sequence {rin,hn)- To show this, choose a convergent subsequence r/^^ of the sequence 
rjn in G+. This is always possible since G+ is compact. Then, by considering the convergent 
sequences Vnigniihm)'^ and one can immediately conclude that hm must be a convergent 
sequence in G^, which proves the claim. To continue, notice from (13. 6p that the effectively 
acting symmetry group is the factor group ((?+ x G'^)/(Z„)diag, where (Z„)diag is formed by the 
pairs {zl2n, zl2n) ^ X G^ with z running over the n^^ roots of unity. We shall demonstrate in 
the proof of Theorem 3.4 that the action of (G+ x G^) / (Z„)diag on P^ is a free action. Moreover, 
Pc is a closed, embedded submanifold of P, as it follows from the definition (I3.10p of Pc and from 
the locally free character of the {G+ x G+)-action on it. Since we have a free and proper action 
on the manifold Pc, the general theory [11] guarantees that Pj-cd — Pc/{.{G^ x G'^)/(Z.„)diag) is 
a smooth symplectic manifold. This is manifest by the model of P^cd constructed below. 

Our goal in what follows is to exhibit a global cross section (a global 'gauge slice') of the 
orbits of x G^ in Pc, which will yield a concrete model of the reduced Hamiltonian systems 
[Py-ed-i^veAiH]^'^). We first present the following lemma, whose proof will also show how to 
construct a convenient global gauge slice. 

Lemma 3.2. The element e'^, q E A in /i2.23\) . andu G C" enter a triple (e^, J, xG^+^{u)) e Pc 
Ii3.10\) if and only if \uj\'^ = 2k, for all j = 1, n and q is regular in the sense of Eg. Ii2.26\) . 

Proof. Let us inspect the moment map constraint for an element of P of the form 

(e^ J, xG^ + ^{u)) with some q e A. (3.14) 
Denoting the projections associated to the decomposition (I2.18P as tt^ and decomposing J as 

J = j+ + J- + j+ + j-^ (3.15) 
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we can spell out the moment map constraint as the conditions 

j+ = -a^c' - 7r+(e(w)), J- = -7r;(e(w)), (3.16) 

and 

7r+(e"^'^'(J)) = (coshadg)(J+ + J+) - (sinhadg)(J; + Jz) = yCZ (3.17) 
Since C" G the vrl projection of equation (13.171) says that 

(coshadg)(J^) - (sinhadg)(J+) = 0, (3.18) 

and its projection requires that 

(coshadg)(J+) - (sinhadg)(Jl) = yC\ (3.19) 

We here used that coshadg maps Ql to Ql and sinhad^ maps to QZl (with — s = =F for 
s = ±). By substituting from fl3.16p into fl3.19p and then taking the scalar product of both 
sides of equation fl3.19|) with an arbitrary T G (12.241) . we obtain the requirement 

(T,e(M)) = VTgTW, (3.20) 

where we also took into account that C' and C" belong to (I2.32p . By using the form of 
M (^TM^ and that of ^(m) ([23SD, we can rewrite (IX^ as the condition 

|Mjf = 2k, Vj = l,...,n. (3.21) 

If fl3.2ip holds, then we can apply the action of the subgroup M^iag of x G+, 

Mdiag := {(/i, /i) G G+ X 6-+ I /i G M} (3.22) 

to replace (without changing q) the element in fl3.14p by an element of the form 

(e'', J, x& + ^{u'')) with the vector u] := V2^, j = 1, n. (3.23) 

We further inspect the moment map constraint for the element (13. 23 p . First looking at the 
block-diagonal components of (I3.19p . we see that the matrix elements {JZ)f.^. are arbitrary real 
numbers for all k = 1, n, and that we must have 

— i/€COsh(gj — Qk) — {Jz) ■ ^, sinh(gj — Qk) = for 1 < j < k < m and m < j < k < n. 

(3.24) 

The last equation can be solved for {J^)j ^ if only if {qj — qk) 7^ for the pertinent indices. 
Next, the block off-diagonal components of (I3.19P can be spelled out as the conditions 

— in cosh(gj + qk) — (JZ) ■ sinh(gj + qk) = for 1 < j < k < m and m < j < k < n, 

(3.25) 

and 

- ix cosh(2gfe) - {JZ)^^^_^f^ sinh(2gfc) = y{C'')k,n+k for k = l, n. (3.26) 

Equation fl3.25p can be solved for {JZ) ■ if and only if {qj + qk) 7^ for the relevant indices. 
Taking into account the assumption (a;^ — y^) 7^ (13. lip and the formula (I2.3ip of C^, equation 
f l3.26p can be solved for ( Jr)^ if and only if 7^ for all k. 

We have seen that equation (13.191) admits a solution if and only if u satisfies (I3.2ip and q is 
regular fl2.26p . The proof is finished by noting that the remaining equation f l3.18p can always 
be solved for if = —tt^{^{u)) is given, since coshad, yields an invertible map on □ 
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Definition 3.3. Suppose that k > and x, y satisfy liS.ll]) . For any q ^ Ac and p E A define 
the function J{q,p) by the formula 

J{q,p) := -x& - i{u^) + L(g,p), (3.27) 

where {u'^)j = \PIk (j = l,...,nj and the matrix elements of L{q,p) = n-{J{q,p)) are the 
following. Firstly, ifl<j<k<m or m < j < k < n, then 

Lj^k = ~Lkj = —Lj^n,k+n = Lk+n,j+n = — i/tcoth(gj — qk) , (3.28) 

Lj^k+n = Lkj+n = —Lj+n,k = —Lk+n,j = — i/t coth(gj + qk). (3.29) 
Secondly, if 1 < j < m and m < k < n, then 

Lj^k = ~Lkj = —Lj^n,k+n = Lk+n,j+n = — i/ttanh(q'j — g^), (3.30) 

Lj^k+n = Lkj+n = ~Lj-f^n,k = ^Lk+n,j = — iKtanh(gj + qk). (3.31) 
Finally, for any 1 < j < m, m < k < n, and 1 < I < n, we have 

Ljj+n = -L,+n,j = -—^^ - coth(2g,), (3.32) 

Lk,k+n = -Lk+n,k = . — T - ix coth.{2qk) . (3.33) 
sinh(2gfc) 

Li^l = -Li+ri,l+n = Pi- (3.34) 



Theorem 3.4. By using the above definition of J{q,p), consider the set 

S = {(e", J{q,p),xC^ + ^{u")) IqeAcpeA}. (3.35) 

The submanifold S G P lies in the constraint surface Pc 113. 10\) and intersects every orbit of 
Gj^ X in Pc precisely in one point. The pull-back Vts of the symplectic form Q 113. 2\) on S is 
given by 

n 

= ^ dpk A dqk. (3.36) 

k=l 

Thus the symplectic manifold {S,VLs) provides a model of the reduced phase space (Pred^^red) 
^3.12\) . which can be identified with the cotangent bundle T*Ac. 

Proof. We know that every g E G can be decomposed according to (12.271) . and Lemma 3.2 
implies that every gauge orbit (i.e. G+ x G^ orbit) in Pc admits a representative of the form 

{e\J,xG^ + i{u^)) with g G A, (3.37) 

where Ac is an open Weyl chamber (for example the one defined in fl2.29p ). Following the 
proof of Lemma 3.2, it is easy to check that J in fl3.37p can be written as J = J{q,p) in (13.271) 
with some p E A. Indeed, the formula (I3.27P was obtained by directly solving the constraint 
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equations listed in the proof of Lemma 3.2. To check that S intersects every gauge orbit only 
once, suppose that we have 

ir]e''h-\7^Jiq,p)7^-\xC'+r]au'')v-') = {e'^\ J{q' ,p'),xC' HW), iv,h) G G+xG+, (3.38) 

for two triples in S. The uniqueness property of the decomposition f l2.27p entails that e'^ = e''\ 
which is equivalent to g = q', and [r], h) = (/i, /i) for some fi G M. Then it follows from the 
second component of the equality in fl3.38p that p = p' holds, i.e., the two representatives of the 
orbit coincide. Incidentally, the equality ^^{u'^)^~^ = C,{u'^) implies that fi E M must belong 
to the centre of M, which is isomorphic to the group Z„ and equals the centre of G. The 
corresponding subgroup (Z„)diag < Miiag acts trivially on P, and hence we can also conclude 
that the factor group (G^ x G^) /{Zn)diag acts freely on the constraint surface Pc- 

We can compute the pull-back of Q fl3.2p on the submanifold S G P, which gives the 
formula (13. 36 p . Since we have seen that S* is a global cross section of the gauge orbits in Pc, it 
follows that {S, fls) represents a model of the reduced phase space (Pred, ^red)- Referring to the 
identification A — A* defined by the scalar product (I2.19p of Q, {S, fls) is symplectomorphic 
to the cotangent bundle T*Ac — Ac x A* equipped with the Darboux symplectic form. 

□ 

Let us recall that a Hamiltonian given by a smooth function on a 2n dimensional symplectic 
manifold is called Liouville integrable if it is contained in a family of n functionally independent, 
globally smooth functions on the phase space whose mutual Poisson brackets vanish and their 
Hamiltonian flows are complete. Now the following result is an immediate consequence of the 
Hamiltonian reduction. 

Corollary 3.5. A family of functionally independent Hamiltonians that are in involution with 
respect to the canonical Darboux Poisson structure on T*Ac is provided by 

= ^tr(J(g,p)2'=), k = l,...,n. (3.39) 

The generalized Sutherland Hamiltonian H{q,p) U.l\) is Liouville integrable, since it obeys 

H{q,p) = \tT{J{q,pf) = Hl^\q,p). (3.40) 



Proof. The reduced Hamiltonians (I3.39P are in involution with respect to the canonical Poisson 
structure derived from Vis (13.360 since the original Hamiltonians Hk (13. 4p are in involution 
with respect to the Poisson structure on (P, f2). By using Definition 3.3, the identity (I3.40p is a 
matter of direct verification. At generic points of the phase space, the Hamiltonians (I3.39P are 
independent, since they start with independent 'leading terms' given by respective homogeneous 
polynomials in pi, The reduction guarantees that the corresponding Hamiltonian flows 

are complete, and thus H]^'^ (and in particular H = ifp'^) is Liouville integrable. □ 
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Finally, let us describe how the flows of the reduced Hamiltonians can be constructed 
from the 'free flows' given in fl3.5p . Take an arbitrary initial value (g(0),p(0)). As a consequence 
of the Hamiltonian reduction, the corresponding solution {q(t),p(t)) of Hamilton's equation for 
Hl'^'^ can be read off from the equality 

(e^W,J(g(t),p(t)),a;C' + eK)) = (3.41) 
= (r/(t)e*^'=e'^(°)/.(t)-\r/(t)J(g(0),p(0))77(t)-\r^(t)(xC' + eK)r7(t)^^) , 

where 

V, = J{q{0),p{0)r-' - i^tr(J(g(0),p(0))2'=-i)l2„ (3.42) 

and {ri(t),h(t)) G G+ x is determined by the condition that the left-hand-side of fl3.4ip 
must belong to the gauge slice S fl3.35p . Thus, finding the solution requires the determination 
of the generalized Cartan decomposition 

etVfceg(o) _ ^(t)-ie^W/i(t), (n{t), q{t), h{t)) eG+xAcX G+, (3.43) 

made unique by the initial condition 7^(0) = h{0) = l2n £ G and continuity in t together with 
the auxiliary condition 

r^{t)au^Ut)-^ = au'^). (3.44) 

Then p{t) obeys 

p{t) = nZ{v{t)J{q{0),p{0)Ut)-') = 7rZ{v{t)L{q{0),p{0)Ut)-'). (3.45) 



If one is interested only in q{t), then a simpler solution algorithm is also available. For this, 
notice that the evaluation of the expression gr{g~^) (with F defined by (12.71) ) for both sides of 
the equality in (13.431) leads to the relation 

e^'/Wn^ = r/(t)e*^'=e25(°)D„e*^'='r/(t)-\ (3.46) 

This means that the entries of the diagonal matrix e^'^^^^D^ are the eigenvalues of the Hermitian 
matrix 

It follows from (13.461) and the forms of Dm (12. 3p and the Weyl alcove (12.291) that the eigenvalues 
of the above Hermitian matrix are all different, and therefore finding q(t) boils down to an 
ordinary diagonalization problem. 

The above algorithm could be particularly useful to analyze the generalized Sutherland 
dynamics, which arises as the k = 1 special case of the reduced systems (I3.39p . The formula 
(11.11) entails that in this case p(t) = q(t). If desired, one could also use the above derivation to 
obtain a Lax representation for the equations of motion, taking J{q,p) in (13.270 . or alternatively 
L{q,p), as the Lax matrix. However, the Lax pair would not give anything substantial to 
our knowledge about the generalized Sutherland model, whose main features follow from its 
realization as a reduction of the free geodesic motion on SU{n,n). 
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4 Conclusion 



In this paper we applied the Hamiltonian reduction approach to a particular many-body system. 
Our derivation of the generalized Sutherland model (11. ip by reduction of the free geodesic 
motion on the group SU{n, n) proves the integrability of the model in the new case of three 
independent coupling constants, and provides a simple algorithm for constructing the solutions. 
This potentially paves the way for future work to analyze the scattering characteristics of the 
model along the lines of the papers [HI US] • The investigation of the quantum mechanics of 
the model (for example by quantum Hamiltonian reduction) is also a challenging problem. 

Another interesting problem is to find duality properties for the generalized Sutherland 
model, which would extend the action-angle dualities of the integrable many-body systems 
studied by Ruijsenaars (see e.g. the review [I5]). This problem exists in general for the Suther- 
land models with two types of particles, whose duality properties are not even known in the An 
case. For the description of dualities in the reduction approach, see also [51 S] and references 
therein. 

Recently [T] new integrable random matrix models have been constructed in association 
with certain integrable many-body systems of Calogero- Sutherland type. It could be feasible to 
extend this correspondence between random matrix models and integrable-many body systems 
to other cases, possibly including generalized Sutherland models with two types of particles. 

We end with a remark on the Lax matrices that can be associated to the model (II. ip . 
Namely, we note that our usage of SU{n, n) as the starting point leads to a 2n x 2n Lax matrix, 
but it should be also possible to derive a (2n + 1) x (2n -|- 1) Lax matrix for the same model, 
with 3 independent couplings, by reduction of the free motion on SU{n + l,n) (cf. Ref. [5]). 
In the case of equal couplings {x = 2y = 2\/2k in (ll.ip ). it is this latter Lax matrix that one 
may expect to obtain directly as well from the standard Lax matrix of the original Sutherland 
model of 2?T, + 1 particles by applying imaginary shifts and restriction to 'mirror symmetric' 
configurations. Although in the reduction approach the role of the Lax matrices is somewhat 
secondary, they are central in other approaches [Ul |T7] . For this reason, we plan to describe 
the alternative Lax matrices of different size and their relationship elsewhere. 
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